Observation of the Charged Hadron Z^±(3900) at ^s = 4170 MeV 
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The BES III Collaboration has reported observation of a resonance, which they call Zc(3900), 
which decays into J /tj} + a charged pion. This charged state was observed as an intermediate res- 
onance in their analysis of e^e~ annihilations into J /xj) + ■k'^tx~ at -/a — 4260 MeV, the peak 
of the vector state Y(4260). We observe this charged state, Z^(3900), at a 6cr significance 
level in the analysis of 586 pb~^ of data taken with the CLEO-c detector at ^(4160). Our 
results, M{Z^) = 3885 ± 5(stat) ± l(syst) MeV, V{Zt) = 34 ± 12(stat) ± 4(syst) MeV, and 
R = cr(e+e" -)■ 7r'FZ^(3900))/cr(e+e" ->• 7r+7r"J/V)) = (35 ± 11)% are in good agreement with 
those reported by BES III. We also present evidence for the observation of the neutral ^^(3900) 
at a 3cr level of significance with parameters in agreement with those of Z^(3900). Because our 
observations are made at the peak of the 2^D\ state of charmonium, we do not find any reason to 
attribute any unconventional properties to the initial state formed by the e^e~ annihilations. 



In a recent posting in the archive [l|, BES III has re- 
ported observation of a charged hadron Zc(3900) with 
mass 3900 MeV, and width ~ 46 MeV, which decays 
into a charged pion and J /ip- This is an important find- 
ing because a charged hadron decaying into a charmo- 
nium state plus a charged meson must contain at least 
four quarks. If confirmed, this observation would her- 
ald the possible existence of a family of other charged 
states in this mass region. Several charged hadrons 
were reported earlier by Belle; in the charmonium region 
Zc(4050, 4250) 0, and ^^(4430) [1, and in the bottomo- 
nium region, Zb(10610, 10650) [J]. So far, these states 
have not been independently confirmed, and many re- 
main controversial. 

The BES observation of Zc was made in e+e^ annihila- 
tion at the vector resonance Y(4260), which is known to 
have unusual characteristics and does not fit in the con- 
ventional spectrum of charmonium states. Clearly, it is 
very important to confirm the existence of .Zc(3900), and 
to also determine if it is populated in e^e~ annihilation 
at a resonance which has well established charmonium 
credentials. We have made a search for Zc in the same 
decay chain as BES 



e^e 



-> 



TT^ZJ 



--J/^, 



but using CLEO-c data taken at -/s — 4170 MeV, on the 
peak of the well-known ipiAlQQ) charmonium resonance. 

We use 586 pb~ of e+e~ collision data taken at 
yi = 4170 MeV at the CESR coUider at the Cornell 
University, with final state particles detected and identi- 
fied in the CLEO-c detector. The CLEO-c detector has 
been described in detail elsewhere [5[. 

We make the event selection for ■0(4160) — )• tt+tt^ J/t/i, 
J/f/) — >■ e^e~ , /x^/i", as in our paper [a] on the search for 
charmonium-like exotic mesons in the radiative decay of 
■0(4160), except that we do not allow any photons. To re- 
construct the TT+TT" J/^0 final state, we select events with 
4 charged particle tracks with zero net charge. Tracks are 
reconstructed in the region with | cos^t^l < 0.93, where 9 



is the polar angle, and are required to be well-measured 
and consistent with originating at the interaction point. 

Charged particle tracks are first identified on the ba- 
sis of their momenta. Leptons (e, /i) from the decay of 
J /^ have momenta > 1 GeV, and pions have momenta 
< 1 GeV, which makes tt/^ separation easy. Pion candi- 
dates are additionally required to have an energy loss in 
the drift chamber {dE/dx) consistent with that expected 
for pions within 3ct. 

Muons are distinguished from the electrons based on 
the variable Ecc/p, where p is the track momentum 
measured in the drift chamber and Ecc is the energy 
deposited in the calorimeter associated with the charged 
particle track. This variable clearly separates electrons, 
which deposit all of their energy in the calorimeter and 
have Ecc/p ~ Ij from muons which pass through the 
calorimeter and deposit only minimum ionizing energy. 
Muons are required to have Eqc/p < 0.25. To re- 
ject backgrounds from photon conversions, we require 
cos(7r+,7r-) < 0.98 and cos(7r±,eT) < 0.98. 

To select fully reconstructed events and improve mass 
resolution, a 4C kinematic fit is performed constrain- 
ing the TT'^TT~J/ip final state to a common vertex with 
Xvcrtcx < 20, and the e~^e~ collision energy and momen- 
tum with xic fit < 20. In the following, we use the mo- 
menta of the charged particles after the kinematic fit. 
To select events containing a J/0 — >■ e"'"e~,/Lt"'"/z~ decay, 
we select events with a e^e^//i+/i^ mass consistent with 
that of the J/tJj within ±20 MeV/ c^. 

The observed J/tp decays are illustrated in Fig. 1. For 
the the well defined peaks of J/tp we obtain for h~^ij,~ 
decay, M{J/tp) = 3096.5±0.5 MeV, N = 137±15 counts 
and a fitted resolution width of 4.2 ± 0.5 MeV, and for 
e+e- decay M{J/iP) = 3096.3 ± 0.5 MeV, TV = 96 ± 
11 counts and a fitted width of 3.9 ± 0.5 MeV. Monte 
Carlo determined efficiencies for decays containing /x^ are 
e = 51%, and for decays containing e^ are e = 43%. We 
estimate the radiative correction factor C = 0.82. The 



cross sections are determined as a{e^e 
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FIG. 1. Invariant mass distributions for J/^ -^ e^e~ (top) 
and J /tjj — > M^M" (bottom) events. The curves show results 
of fits with gaussian peaks and constant background. The 
PDG average M{J/i)) = 3096.916 ± 0.011 MeV [1]. 



N/C e C. They are cr = 9.5 ± 1.0 pb (^i^), and tr = 7.9 ± 
0.9 pb (e^), with the average value of 8.6 ± 0.7(stat) pb. 
This is nearly 1/7 th of the cross section observed by BES 
in decays from Y(4260). 

We construct distributions of the maximum of 
M{iT+J/ip) and M{n-J/^), labeled Mmax(7r=^^/V'), the 
same observable as BES, in two regions defined as the 
J/i; resonance region, M(e+e~, ^+/^") = 3097±20 MeV, 
and the sideband regions, M{e'^e—,iJ.~^fj.~) = 3030 ± 
30 MeV and M(e+e-, /x+At") = 3170 ± 30 MeV. We find 
that there is no peaking in the distribution for the side- 
bands, and that its shape is identical within statistical 
uncertainties with the shape of the phase space distri- 
bution. Fig. 2 shows the Mi„ax(7i'^ J/i/') distribuion for 
the resonance region. The enhancement at ^ 3880 MeV, 
which we attribute to the Zc resonance, is clearly visible. 



FIG. 2. Distribution of A/max (t J/ ■0), as observed in the 
decays (top) ^(4160) — > -k^tt' J/ip and (bottom) ^(4160) — > 
TT TT J/ip- The points with errors bars are data. The his- 
togram shows the results of the fit with the hashed part show- 
ing the Breit-Wigner function convolved with a gaussian res- 
olution function for the Zc resonance, and the shaded part 
showing the MC-determined shape for phase-space J/tpnn de- 
cays. 



The distribution was fitted assuming the phase space MC 
shape for the background and a peak whose shape was 
obtained by convolving the MC determined instrumen- 
tal width of 3.8 MeV with a Breit-Wigner resonance. 
The normalization of the background, the position of 
the peak, its width, and its magnitude were kept free 
in the fit. The result of the fit was number of counts, 
iV = 81 ± 20, mass M{Zc) = 3885 ± 5(stat) MeV, and 
width T{Zc) = 34 ± 12(stat) MeV. The x^/d-o.f. of the 
fit is 35/37. The likelihood-determined significance of the 
peak is 6.1(7. 

We evaluate systematic uncertainties as follows. Our 
fitted J/il} masses differ from the known M{J/ip) = 



TABLE I. Summary of fit results and comparison with results from BES III ^. 



Vs C N(J/i) -^ a{e+e- -^ N{Z^) signif. M{Z^) V{Z^) R 

GeV pb-^ ^+/i-/e+e-) ttttJ/t/;), pb g MeV MeV % 



BES III TT+TT" 4.26 525 882 ± 33/595 ± 28 62.9 ±1.9 307 ± 48 >8 3899 ± 6 

Present tt+tt" 4.17 586 137 ±15/ 96 ±11 8.6 ± 0.7 81 ± 20 6 3885 ± 5 

7rV° 4.17 586 40 ± 8/ 29 ± 5 5.7 ± 0.8 17 ±10 3 3907 ±12 



46 ± 22 22 ± 3 
34 ±13 35 ±11 
34 ±29 25 ± 15 



3096.916 ± 0.011 MeV by less than 0.5 MeV, which im- 
plies that the systematic uncertainty in our mass cal- 
ibration is less than 0.5 MeV. To estimate the uncer- 
tainties due to the choice of the phase space shape for 
the background, we consider the background shape used 
by BES, a/{x - 3.6)'' + c + dx [I]. An equally good fit 
is obtained with y^ /d.o.f. = 33/34. This leads to sys- 
tematic uncertainties of ±0.2 MeV in mass, ±4 MeV 
in width, and ±17% in cross section. We estimate the 
uncertainty in the MC-determined resolution width by 
varying it by ±10%; it leads to systematic uncertain- 
ties of < ±0.1 MeV in mass, ±0.5 MeV in width, and 
< 1% in cross section. We estimate the systematic 
uncertainties due to track reconstruction to be 1% per 
track. We estimate a 1% uncertainty due to the lu- 
minosity measurement [7|. The branching fractions for 
J/ip — >■ e+e~,/x+/.t~ have an uncertainty of 1% [8|. In 
quadrature, all these add to a systematic uncertainty of 
±1 MeV in mass, ±4 MeV in width, and 18% in the cross 
section. 

As Table I shows, our final results, including sys- 
tematic uncertainties are M{Zc) — 3885 ± 5 MeV and 
r(Zc) = 34± 13 MeV. These are in good agreement with 
those reported by BES^ M{Zc) = 3899 ± 6 MeV, and 
T{Zc) = 46 ± 22 MeV fl. We obtain the ratio 



R = a{e+e- -^ n^ Z^{3900j)/cr{e+e- 



-^ vr+TT J/ip) 
(35 ±11)%, 



which is also in good agreement with R = (21.5 ± 3.3)% 
reported by BES [1] . 

The Z^ state decaying into J/ipn^ has isospin / — 
1. If the pion is in the S-state, with respect to J/'ip, 
J^'-^ {Zc) = 1"' , which is not allowed for a qq meson. As 
has been noted by Maiani et al. [9|, it is interesting to look 
for the neutral member of this "exotic" isospin triplet 
in the decay to tt^tt" J/-)/). Taking account of isospin 
and the nearly factor two smaller efficiency for detect- 
ing 7r*'7r'^ compared to tt+tt^, we expect to observe only 
~ 60 J/i/iTT^TT^ events, and ^ 20 Z°. Despite the limited 
expected statistics, we have made this measurement. The 
observed distribution of J /ipTi'^ events is shown in Fig. 3. 
The distribution is fitted in exactly the same manner 
as in Fig. 2. The results are N{Z°) = 17 ± 10 counts. 



Af (ZO) = 3907 ± 12 MeV, r(ZO) = 34 ± 29 MeV, and 
R — (25 ± 15)%. The likelihood determined significance 
of the Z° peak is 3(t. These results are also listed in 
Table 1, and are in agreement with those for Z^ within 
their statistical errors. 

To summarize, our observations provide independent 
confirmation of the existence of the Z^(3900) state, and 
provide new evidence for the existence of the neutral 
member Z°(3900) of this isospin triplet. It is important 
to note that our observations of Z^''^(3900) are based 
on e^e" annihilations at -y/s — 4170 MeV, on the peak 
of the well-established 2^Di vector state of charmonium, 
and therefore we do not need to attribute any unconven- 
tional properties to the initial state, as may be invoked 
for the observation of Zf at Y(4260). 

This investigation was done using CLEO-c data, and 
as members of the former CLEO Collaboration we thank 
it for this privilege. This research was supported by the 
U.S. Department of Energy. 

Note added — As we were preparing to submit this 
paper, we became aware of a paper from the Belle 
Collaboration that also reports on the Z^(3900) from 
Y(4260) [i3|. 
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